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ABSTRACT 



Aims. The blue compact dwarf galaxy IIZw71 is catalogued as a probable polar-ring galaxy, and along its long axis it has several very 
luminous knots showing recent episodes of star formation. Our main aim is to study the physical properties, the stellar content, and 
, the kinematics in the brightest knots of the polar ring 

Methods. We carried out long-slit spectroscopic observations of the polar ring in the spectral range 3500 - 10000 A taken with the 
William Herschel Telescope (WHT). The spectroscopic observations complemented by the available photometry of the galaxy in 
narrow Ho- filters. 

Results. We measured the rotation curve of the ring, from which we infer a ratio M/L B w 3.9 inside the star-forming ring. We 
measured the auroral [Om] line in the two brightest knots, allowing us to measure oxygen, sulphur, nitrogen, argon, and neon chemical 
abundances following the direct method. Different empirical calibrators were used to estimate the oxygen abundance in the two faintest 
knots where the temperature sensitive lines could not be measured. The metallicities obtained are very similar for all the knots, but 
£S) ■ lower than previously reported in the literature from integrated spectra. The N/O abundance, as derived from the N2O2 parameter 

(the ratio of the [Nn] and [On] intensities), is remarkably constant over the ring, indicating that local polution processes are not 
■^j- ' conspicuous. Using synthetic stellar populations (SSPs) calculated with the code STARLIGHT, we studied the age distribution of the 

\ stellar populations in each knot, finding that in all of them there is a combination of a very young population with less than 10 Myr, 

, responsible for the ionisation of the gas, with other populations older than 100 Myr, probably responsible for the chemical evolution 

■ of the knots. The small differences in metallicity and the age distributions among the different knots are indicative of a common 

chemical evolution, probably related to the process of interaction with the companion galaxy IIZw70. Finally, we calculated the SFR 
in the different knots from the Ha luminosities. The combined SFR rate for the ring amounts to half of the integrated SFR for this 
galaxy reported by previous authors. The average surface SFR is also higher than the average reported value of HII in polar rings by 
more than an order of magnitude. 
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C3 ' 1 . Introduction tant to study the cases taking place in the nearby Universe now 

and the links between dwarf interactions and the star formation 

The intense processes of star formation occurring in the lumi- history of BCDs 
nous knots of HII galaxies make them easily detected in surveys 

based on strong emission lines. In fact, these knots present spec- Although IIZw71 is catalogued as a BCD galaxy it is, in 
tra similar to those of HII regions, so it is possible to use them to fact) characterised by several very luminous Ha knots distributed 
ascertain the basic properties of the host galaxy, such as extinc- a i ong a ring that is rotating around the host galaxy. It has been 
tion, metallicity, star formation rates, and from them derive the catalogued as a probable polar-ring galaxy (B17) in the Polar- 
chemical evolution and the star formation history of the whole Ring Catalogue (Whitmore et al., 1990). Polar ring galaxies 
galaxy. (PRGs) are systems with two kinematically separated compo- 
HII galaxies are a subset of blue compact dwarf galaxies nents. The central component (the host galaxy) is usually a 
(BCDs), which are selected by their compact aspect and blue lenticular galaxy or occasionally an elliptical galaxy. The other 
colour. This type of object is thought to be more frequent in a component, the polar ring, follows a approximately circumpolar 
younger Universe and they are possible building blocks for the orbit around the host and it is characterised by the presence of 
largest galaxies that we can detect at low redshifts (Kauffmann sta rs, molecular gas ,and dust, inside another bigger ring com- 
et al., 1993). If interactions among dwarf irregular galaxies was poS ed of neutral hydrogen. Thus, this ring becomes an appro- 
a basic mechanism of galaxy formation in the past, it is impor- pr iate place for star formation. It is thought that these objects 

are formed as a consequence of the interaction between galax- 

Send offprint requests to: E. Perez-Montero ies with a small impact parameter (Bournaud & Combes, 2003). 
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In the case of IIZw71, there are proofs from interferometric ob- 
servations of interaction with IIZw70, another BCD. In fact, al- 
though clearly separated on optical images, both galaxies share a 
common HI envelope with a gaseous bridge or streamer connect- 
ing both structures (see Figure 3 of Cox et al., 2001). This points 
to an ongoing interaction between the two galaxies. Therefore 
this system becomes an ideal scenario for studying the effects of 
interactions in the formation and evolution of BCDs in the Local 
Universe. A distance of 18.1 Mpc to the system is adopted by 
Cox et al. (2001), taking a value of H = 75 Km s~' Mpc" 1 . This 
implies a linear scale of 90 pc/arcsec on the sky. 

In a single long-slit exposition, we have observed the main 
bursts of star formation along the direction of the polar ring of 
IIZw71. The analysis of the obtained spectra and the available 
photometry in several bands were used to make a comparative 
analysis of the different bursts of star formation. In the next sec- 
tion we describe the long-slit spectroscopic observations of the 
polar ring of IIZw71. In sect. 3, we present the results of our 
study, including the determination of physical properties, such as 
electron density and reddening in the four observed knots, and 
of chemical abundances in the two brightest knots. In sect. 4, 
we discuss our results including a study of the kinematics of the 
polar ring, the determination of metallicity in both the faintest 
and the brightest knots, the reddening, the stellar properties by 
means of fitting synthesis stellar populations, and the star forma- 
tion rates in the individual knots using the available Ha photom- 
etry of this galaxy. Finally, our conclusions are presented in the 
last section. 

2. Observations and reduction 

The long-slit spectrophotometric observations of IIZw71 were 
obtained using the ISIS double-beam spectrograph mounted 
on the 4.2m William Herschel Telescope (WHT) of the 
Isaac Newton Group (ING) at the Roque de los Muchachos 
Observatory on the Spanish island of La Palma. They were ac- 
quired on 2005 July 8 during one single night's observing run 
and under photometric conditions, with an average seeing of 
0.7 arcsec. The EEV12 and Marconi2 detectors were attached 
to the blue and red arms of the spectrograph, respectively. The 
R600B grating was used in the blue covering the wavelength 
range 3670-5070 A (centred at /1 C =4370A), giving a spectral 
dispersion of 0.45ApixeU'. On the red arm, the R316R grat- 
ing was mounted in two different central wavelengths provid- 
ing a spectral range from 5500 to 7800 A (A c = 6650 A) and 
from 7600 to 9900 A (A c = 8750 A) with a spectral dispersion 
of 0.86ApixeL'. To reduce the readout noise of our images, 
the observations were taken with the 'SLOW CCD speed. The 
pixel size for this set-up configuration is 0.2 arcsec for both spec- 
tral ranges. The slit width was ~1 arcsec, which, combined with 
the spectral dispersions, yields spectral resolutions of about 1 .0 
and 3.5 A FWHM in the blue and red arms, respectively. The 
instrumental configuration, summarized in TableQ] was planned 
to cover the whole spectrum from about 3500 to 10000 A , at 
the same time providing a moderate spectral resolution. This 
guarantees the simultaneous measurement of the nebular lines 
of [On] AA 3727, 3729, and [Sin] AA 9069,9532 A at both ends of 
the spectrum, in the very same region of the galaxy. 

To cover all knots along the polar ring at the same expo- 
sure, all the spectra were taken at PA= 30°, an average of 56° 
off paralactic angle, and the spectra are thus somehow affected 
by atmospheric differential refraction. Taking the mean air mass 
into account during the observations of this object (« 1.25) and 



Table 1. WHT instrumental configuration 



Spectral range Disp. FWHM Spatial res. Exposure Time 

(A) (Apx- 1 ) (A) ("px-') s 

3670-5070 0.45 1.0 0.2 4x 1200 

5500-7800 0.86 3.5 0.2 2x 1200 

7600-9900 0.86 3.5 0.2 2x 1200 



following the curves by Filippenko (1982), we calculated that 
the displacement between [On] 3727 A and [Sin] 9069 A, which 
are the emission lines the most separated in wavelength, is no 
greater than 0.9 arcsecs. 

Several bias and sky flat field frames were taken at the be- 
ginning and at the end of the night in both arms. In addition, 
two lamp flat fields and one calibration lamp exposure were per- 
formed at each telescope position. The calibration lamp used 
was CuNe+CuAr. The images were processed and analysed with 
IRAlfl routines in the usual manner. This procedure includes 
the removal of cosmic rays, bias subtraction, division by a nor- 
malised flat field, and wavelength calibration. Typical wave- 
length fits were performed using 30-35 lines in the blue and 
20-25 lines in the red and polynomials of second to third or- 
der. These fits were made at 117 different locations along the 
slit in each arm (beam size of 10 pixels) obtaining rms residuals 
between ~0.1 and ~0.2pix. In the last step, the spectra were cor- 
rected for atmospheric extinction and flux-calibrated. For both 
arms, two standard star observations were used, allowing a good 
spectrophotometric calibration with an estimated accuracy of 
about 5%. 

The left panel of Figure Q] shows the spatial distribution of 
the Ha flux along the slit. The emission line profile in this plot 
has been generated by collapsing 1 1 pixels of the spectra in the 
direction of the resolution at the central position of the line in the 
rest frame, A 6563 A. A continuum profile was also generated by 
collapsing 1 1 resolution pixels centred at A 6593 A . The differ- 
ence between the two, corresponds to the pure Ha emission. It 
can be see from the profiles that in knot C, the most intense one, 
the peaks of the continuum and line emission do not coincide 
spatially but are separated by 1.6 arcsec. Also the width of the 
continuum profile is larger than that of the line emission. This 
suggests that most of this continuum is coming from the host 
galaxy behind the knot. 

The regions of the frames to be extracted into one- 
dimensional spectra corresponding to each of the knots were 
selected on the Ha line and are marked by horizontal lines in 
the left panel of Figure Q] In the right panel of the same figure 
we show the Ha image and continuum contours in the R-band 
from Gil de Paz et al. (2003) in order to illustrate the position of 
the bursts of star formation and the relative position of the host 
galaxy. We also show elliptical regions in Ha taken to measure 
the total Ha flux of the knots for an aperture correction of the 
slit measurements. The different emission knots are labelled in 
both panels. 

Figure [2] shows the blue and red parts of the spectra of the 
four observed knots (labelled from A to D), which have been 
extracted from the slit. Knot C, as expected, shows the greatest 
contribution from an underlying stellar population. 



1 IRAF: the Image Reduction and Analysis Facility is distributed by 
the National Optical Astronomy Observatories, which is operated by 
the Association of Universities for Research in Astronomy, In. (AURA) 
under cooperative agreement with the National Science Foundation 
(NSF). 




Fig. 1. In left panel, spatial profile of the light distribution along the slit for the observed Ha emission (dashed line), the adjacent 
continuum (thick green line), and the difference between the two (thin blue line). The horizontal lines show the spatial pixels that 
have been compressed to produce the one dimensional spectra corresponding to each emission knot. In the right panel, we show 
the Ha image from Gil de Paz et al. (2003) with the identification of the observed knots and the regions for which the broad band 
magnitudes and the narrow Ha flux have been measured. In the same image, we show the position of the slit for the observations 
described in the text and the elliptical regions measured to correct for aperture effects in each knot. We show the R-band contours 
as well, as taken from the same source, which show the position of the host galaxy. 



3. Results and discussion 

3.1. Line intensities and reddening correction 

The presence of a conspicuous underlying stellar population de- 
presses the Balmer emission lines and does not allow the mea- 
surement of their fluxes with an acceptable accuracy (Diaz, 
1988). To account for this, we subtracted from the observed 
spectra the best-fitting found by the spectral synthesis code 
STARLIGHT (Cid Fernandes et al. 2004, 2005; Mateus et al. 
2005). STARLIGHT fits an observed continuum spectral energy 
distribution using a combination of multiple simple stellar popu- 
lation (SSPs; also known as instantaneous burst) synthetic spec- 
tra using a x 2 minimisation procedure. We chose for our anal- 
ysis the SSP spectra from Bruzual & Chariot (2003), based on 
the STELIB library of Le Borgne et al. (2003), Padova (1994) 
evolutionary tracks, and a Chabrier (2003) initial mass function 
between 0.1 and 100 M . We fixed the metallicity to Z = 0.004 
(= 1/5 Z ) for all models, which is the metallicity in the Bruzual 
& Chariot (2003) library closest to those corresponding to the 
oxygen abundances measured in knots B and C, as explained in 
Section 3.3. The code simultaneously finds the ages and the rel- 
ative contributions of the different present SSPs and the average 
reddening. The reddening law from Cardelli, Clayton & Mathis 
(1989) was used. Prior to the fitting procedure, the spectra were 
shifted to the rest frame, and sampled again to a resolution of 1 
A. Bad pixels and emission lines were excluded from the final 
fits. The percentage of the mean deviation over all fitted pixels 
between observed and model spectra ranges between 7.3% in 



2 The STARLIGHT project is supported by the Brazilian agencies 
CNPq, CAPES and FAPESP and by the France-Brazil CAPES/Cofecub 
program 



Table 2. Equivalent widths of the hydrogen recombination lines 
H/3, Hy and H5 as measured directly in the spectra of the ob- 
served knots once the underlying stellar populations have been 
removed, with the corresponding correction factors. 



Knot 


A 


B 


C 


D 


EW(Hq-) (A) 


270+30 


130+10 


40±2 


140±5 


f c 


1.00 


1.00 


1.05 


1.00 


EW(Ha) t . (A) 


270+30 


130+10 


42±2 


140±5 


EW(Hfi) (A) 


34.9±4.5 


30.0±1.7 


10.7+1.1 


32.2+4.7 


fc 


1.01 


1.03 


1.13 


1.00 


EW(H/?) C (A) 


35.2±4.7 


30.9+1.8 


12.1+1.3 


32.2+4.9 


EW(Hy) (A) 


17.2+1.7 


12.8±0.9 


5.0±0.9 


10.8+1.7 




1.03 


1.03 


1.19 


1.00 


EW(Hy) c (A) 


17.7+1.3 


13.2+1.0 


6.0+1.0 


10.8+0.9 


EW(H<S) (A) 


8.1±0.6 


6.2±0.7 


1.8+1.0 


8.0+0.8 


fc 


1.09 


1.13 


1.24 


1.03 


EW(Hd) c (A) 


8.8+0.6 


7.0±1.0 


2.3±1.0 


8.2+1.0 



knots A and D, which have the lowest signal, to 2.9 % in knot B 
and 1.9% in knot C. 

Figure|3]illustrates the results of the model fitting to the spec- 
trum of knot C. The upper panel shows the observed spectrum 
and the model spectrum. The lower panel shows the subtraction 
of the fitted stellar continuum. We measured the Balmer emis- 
sion line fluxes over a linear continuum derived by hand from 
the observed spectra. We then compared these with the same 
measurements over the continuum derived from the model SSPs. 
This allows the equivalent widths to be corrected by the pres- 
ence of absorption stellar features in these lines that cause the 
measured EWs to be underestimated. In Table [2] we summarise 
our results for the four brightest Balmer emission lines. For each 
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Fig. 2. Blue and red spectra up to 7200 A for each of the extracted knots of IIZw7 1 . 
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Fig. 3. Results of the STARLIGHT fit to the continuum spectral 
energy distribution for knot C, the one most affected by under- 
lying stellar population contribution. The lower panel shows the 
emission line spectrum once this contribution has been removed. 



Balmer line, the first row gives the measured equivalent width 
of the emission line, while the third row gives the corresponding 
value after the removal of the underlying absorption. The second 
row gives the ratio between the two. Knot C is the most affected 
indicating a larger contribution by the underlying stellar popu- 
lation, although the correction for the Hy8 line only amounts to 
13%, reaching 24% for the HS line. For the rest of the knots, only 
B presents a significant contribution of 13 % for the H5 line. 



The emission line fluxes were measured on the spectra us- 
ing the splot task of IRAF and are listed for the four observed 
knots in Table [3] Column 1 of the table shows the wavelength 
and identification of the measured lines. The observed emission 
line fluxes, F(A) (in units of H/3 flux = 1000) with their corre- 
sponding errors, are given in columns 2, 4, 6, and 8 for Knots A, 
B, C, and D respectively. Only the fluxes of the Balmer recom- 
bination lines were measured on the subtracted spectra, while 
the collisional emission lines were measured on the uncorrected 
spectra in order to reduce the uncertainties due to the subtraction 
of the continuum. 

We used two different ways to integrate the flux of a given 
line. (1) In the case of an isolated line, the intensity was calcu- 
lated by integrating between two points given by the position of 
a local continuum placed by eye. (2) In the case of the two lines 
of [On] at AA 3727, 29 A, which are partially blended, a multiple 
Gaussian fit procedure to estimate individual fluxes was used. 
This procedure is illustrated in Figure [5] The statistical errors 
associated with the measured emission line fluxes were calcu- 
lated using the expression 

a-, = cr c 7V 1/2 [l +EW/(NA)] 1/2 

where 07 is the error in the measured line flux, <r c represents 
the standard deviation in a box near the measured emission line 
and stands for the error in the continuum placement, tV is the 
number of pixels used in the measurement of the line flux, EW 
is the line equivalent width, and A is the wavelength dispersion 
in Angstroms per pixel. 

The reddening coefficients (c(Hy6) were calculated from the 
measured Balmer decrements, F(A)/F(H/3), adopting the galac- 
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Table 3. Relative observed and reddening corrected line intensities [F(H/J)=I(Hj8)=1000] for the different observed knots 
of IIZw71. For each knot we also give the flux of H/3, once reddening corrected and after subtraction of the underlying 
absorption and the constant of reddening, with their corresponding errors. 



A (A) A B C D 





F(X) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


3726 [Oii]^ 


2018+62 


2070+240 


1545+43 


1667+132 


1497±60 


1561+133 


1225±107 


1552±227 


3729 [Onf 


2936±289 


2767±377 


2271 ±61 


2450±163 


1990±74 


2076±156 


1758+150 


2227±290 


3869 [Nem] 


321 ±37 


328±93 


242±23 


258±50 


297±23 


307±53 


183+25 


224±68 


4102 m 


303±29 


307±75 


242+15 


254±40 


257±19 


263 ±42 


233±32 


270±67 


4340 Hy 


584+91 


570+117 


466±13 


480±43 


483+31 


491 ±52 


417±43 


458±75 


4363 [Ora] 






29±5 


30+11 


43±10 


44±16 






4861 H/3 


1000±39 


1000±39 


1000+37 


1000±37 


1000±47 


1000±47 


1000±118 


1000±118 


4959 [Om] 


734±42 


733±54 


818+27 


814±34 


937±34 


934±41 


850±75 


837±80 


5007 [Om] 


2101±69 


2096±99 


2411 ±67 


2393±79 


2770±95 


2759±106 


2608±119 


2550±223 


6300 [Oi] 


83±19 


81±30 


92±7 


87+16 


77±6 


74±15 


100+13 


84±25 


6548 [Nil] 


101±19 


99±19 


92±7 


86±8 


93±13 


90±17 


100+13 


82+15 


6563 Ha 


2752+126 


2693+269 


2895±93 


2713±161 


3000±120 


2894±183 


3333±324 


2724±320 


6584 [Nil] 


303±21 


296±24 


250±18 


234±21 


230±14 


222±18 


250±15 


204±19 


6717 [Sn] 


514+41 


502±59 


405±17 


378±29 


297±12 


285±24 


267±15 


215±29 


6731 [Sn] 


339±30 


332±47 


292±12 


272±24 


217+14 


208±24 


200±9 


161±24 


7137 [Arm] 


138±19 


134±46 


68±7 


63+18 


77±6 


73±20 






9069 [Sin] 


92±28 


88±42 


124+15 


1 10+43 


107±14 


100±43 


250±43 


174±80 


F(H/?)(10- 1:, ergseg- 


1 cm- 2 ) 0.89 


+0.05 


3.61±0.10 


4.19±0.10 


3.49±0.10 


c(H/J) 


<0.05 


0.09z 


=0.03 


0.05= 


t0.02 


0.28±0.02 



° Partially blended (see text and Figure[5]l 



tic extinction law of Cardelli et al. (1989) with R v -3.l. A 
least square fit of the measured decrements to the theoretical 
ones, (F(A)/F(H/3))o, computed based on the data by Storey 
& Hummer (1995), was performed that provides the value of 
c(H/3). The theoretical Balmer decrements depend on electron 
temperature and density. We used an iterative method to esti- 
mate them, taking as starting values those derived from the mea- 
sured [Sn] AA 6717,6731 A and [Ora] AA 4363, 4959, 5007 A 
line fluxes. 

In Figure |4] we show the actual fits for each of the observed 
knots. The Galactic extinction for this object is negligible, so we 
can be sure that internal reddening is the unique source of ex- 
tinction. All the observed knots present positive values of c(H/T), 
except knot A in which the extinction is consistent with zero 
value. The highest value of c(H/3) is found for knot D. 

The values of c(H/3)) and their corresponding errors, consid- 
ered as the uncertainties of the least square fits, are listed at the 
bottom of Table [3] for the observed knots, together with the red- 
dening corrected H/3 intensity. The corrected intensities for the 
observed lines relative to H/3, 1(A), are given in columns 3, 5, 7, 
and 9 of the same table. In all cases, the reddening correction 
was done relative to the closest Balmer recombination line. 

3.2. Physical conditions of the gas 

The physical conditions of the ionized gas, including electron 
temperatures and electron densities, were derived from the emis- 
sion line data using the same procedures as in Perez-Montero & 
Diaz (2003), based on the five-level statistical equilibrium atom 
approximation in the task temden of the software package IRAF 
(De Robertis, Dufour & Hunt, 1987; Shaw & Dufour, 1995). 
The atomic coefficients used here are the same as in Hagele et 
al. (2006). We took as sources of error the uncertainties associ- 
ated with the measurement of the emission-line fluxes and the 
reddening correction and we propagated them through our cal- 




-0.3 -0,2 -0,1 0,0 0,1 0,2 -0,3 -0,2 -0,1 0,0 0,1 0,2 

f(A.) - f(HP) 

Fig. 4. Representation for each observed knot of the ratio of the 
expected and observed ratios of the Balmer hydrogen lines and 
the extinction curve. The slope of the obtained fit is the constant 
of reddening. 



culations. Electron densities were derived in the four knots from 
both [On] AA 3726 / 3729 A and [Sn] AA 6717/6731 A line ratios, 
which are representative of the low-excitation zone of the ion- 
ized gas. The values for each of the knots are listed in Table [4] 
In all cases they provide upper limits, which are lower when the 
[On] lines are used. In the case of these lines, however, the spec- 
tral dispersion of our data do not allow total resolution of the 
lines, and they were deconvolved by fitting two Gaussian com- 
ponents as shown in Figure [5] The upper limits for the electron 
density are lower in all cases than the critical value. 

For knots B and C, for which the intensity of the [Om]/l 
4363 A was measured with sufficient signal-to-noise (see Figure 
[6] for a more detailed plot of this line in these knots), the 
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Fig. 5. Deblending of the [On] AA 3727, 3729 A lines. Their ratio 
allows to put lower limits to the electron density in the observed 
ionized regions. 



Table 5. Electron temperatures, ionic and total abundances 
for the knots B and C of IIZw7 1 . 



Knot 


B 


C 


T([Om]) 


Iz400±l /00 K 


13 /00+zOOO K 


T([Sm]t 


izyuu±i.3uu i\ 


IZoUU+ljUU IX 


11600+1100 K 


13200+1200 K 


tZ+log^C I ri ) 


1 .0 / +U.ZJ 




lZ+log(0 1 H ) 


/. 0.3+0.19 


H CI i A 1 ft 

7.57+0.19 


lZ+10g(U/Hj 


t o^_i_n o 1 

/.yj±u.zt 


n on_i_n o 1 

/ .yu+u.zt 


IZ+lOg^O j tl ) 


c o/i_t_n 1 1 


J.OJxU.l j 


lZ+log(o 1 tl ) 


a nn_i_n oa 


£ 55 j_n OA 

j.oj+u.zo 


T/"^T7/ V + i o2+\ 

K,r(i +0 ) 


1 IC ifl 11 
1. 15+U. 13 


i 1 1 i n n 
1. 14+U. 1Z 


12+log(S/H) 


6.33+0.24 


6.20+0.24 


12+log(/V + //f + ) 
log(N/0) 


5.87+0.28 
-1.48 ±0.37 


5.63+0.23 
-1.68 ±0.32 


12+log(/Ve- + //f + ) 


7.12+0.29 


7.06+0.28 


ICF(2Ve 2+ ) 


1.22+0.13 


1.23±0.14 


12+log(Ne/H) 


7.19+0.35 


7.13±0.35 


12+\og{Ar 1+ /H + ) 


5.61+0.22 


5.57±0.19 


ICF(Ar 2+ ) 


1.16+0.03 


1.16±0.04 


12+log(Ar/H) 


5.65+0.26 


5.60±0.26 



" From a relation with T([Om]) based on photoionization models 
b From an empirical relation with T([Om]) 



5e-16 



IIZw71 - Knot B 



[OIII1 

4363 A 





IIZw71 - Knot C 




I 4320 4340 4360 4380 4400 4420 
Wavelength (A) 

Fig. 6. Detail of the rest-frame spectra of the knots B and C 
around the [Om] 4363 A emission line. 

Table 4. Electron densities in particles per cm 3 of [Sn] and [On] 
measured for each of the knots. 



Knot 


n([On]) 


n([Sn]) 


A 


<200 


<230 


B 


<85 


<160 


C 


<180 


<210 


D 


<190 


<350 



[Om] electron temperature was derived directly from the ratio 
(I(4959A)+I(5007A))/I(4363A). Then we derived T([On]) from 
T([Om]) using the relations based on the photoionization models 
described in Perez-Montero & Diaz (2003), which take explicitly 
nto account the dependence of T([On]) on the electron density, 
and T([m]) has been estimated from the empirical relation: 

r([Sm]) = (1.19 ± 0.08) T([Om\) - (0.32 + 0.10) 



found by Hagele et al. (2006) for HII galaxies. The values of the 
derived electron temperatures for these two knots are listed in 
Tabled 



3.3. Ionic and total chemical abundances in knots B and C. 

The oxygen ionic abundance ratios, + /H + and 2+ /H + , were 
derived from the [On] AA 3727,3729 A and [Om] AA4959, 
5007 A lines, respectively using the appropriate electron temper- 
atures for each ion. The total abundance of oxygen was derived 
assuming 

o ^ + + 2+ 
H ~ IF 

We derived S + abundances from the fluxes of the [Sn] emis- 
sion lines at XI 6717 ', 6731 A assuming that T([Sn]) * T([On]) 
and S 2+ abundances have been derived from the fluxes of the 
near-IR [Sm] A 9069 line and the estimated value of T([Sm]). The 
total sulphur abundance was calculated using an ionization cor- 
rection factor (ICF) for S + +S 2+ according to Barker's (1980) for- 
mula, which is based on the photoionization models of Stasihska 
(1978): 



ICF(S + +S 2+ ) 



O 



2+ 



+ + 2+ 



-1/a 



where a = 2.5 provides the best fit to the scarce observational 
data on S 3+ abundances (Perez-Montero et al. 2006). Taking this 
ICF as a function of the ratio 2+ /0 instead of + /0 reduces the 
propagated error for this quantity. The ICF found are similar for 
the two knots and increase the derived S + +S 2+ by about 15 % . 

The ionic abundance of nitrogen, N + /H + , was derived from 
the intensities of the /1/16548, 6584 A lines assuming that 
T([Nn]) * T([On]). Then, the total N abundance was derived 
under the following assumption: 



N 
O 



o + 
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Neon is only visible in the spectra by means of the [Nem] 
emission line at /13869 A. For this ion we took the electron tem- 
perature of [Oni], as representative of the high excitation zone. 
The total abundance of neon was calculated using the following 
expression for the ICF (Perez-Montero et al., 2007): 

7+ 2+ + + 2+ 

ICF(Ne^) * 0.753 + 0.142 • + 0.171 • 

This formula considers the overestimate of Ne/H in objects with 
low excitation, where the charge transfer between 2+ and H° 
becomes important (Izotov et al., 2004). 

The only accesible emission lines of argon in the optical 
spectra of ionized regions correspond to Ar 2+ and Ar 3+ . In 
the spectra of knot B and C, only [Arm] A 7136 A was mea- 
sured and the abundance of Ar 2+ was calculated assuming that 
T([Arm]) «T([Sm]) (Garnett, 1992). The total abundance of Ar 
was then calculated using the ICF(Ar 2+ ) derived from photo- 
ionization models by Perez-Montero et al. (2007): 

ICF(Ar^) = 0.749 + 0.507 ■ (l - ^^t)+ 

/ o 2+ r 1 

+0.0604- 1 -— — — ~— I 

\ + + 2+ / 

The calculated ICF for Ar and Ne are the same in the two knots 
within the errors. The ionic and total abundances for the ob- 
served species in knots B and C are given in Table |5J along with 
their corresponding errors. 

3.4. Metal content of the polar ring 

The emission-line spectra of the four star-forming knots in 
IIZw71 are remarkably similar, implying similar values for ion- 
ization parameter, ionization temperature, and chemical abun- 
dances. We derived the ionization parameters from the ratio of 
the [On] and [Ora] lines according to the expression given in 
Diaz et al. (2000). They are similar in all the knots ranging from 
6.42xl0~ 4 for knot A to 9.67xl0~ 4 for knot C. Using these val- 
ues, the corrected Ha fluxes, and the sizes of the regions from 
Ha images we can calculate the density of the emitting gas (Diaz 
et al. 1991). This is similar for the four knots with a value of 
about 20 particles • cm 4 , consistent with the upper limits we 
derived from the ratio of the [Sn] lines and providing filling fac- 
tor for the gas of a few times 10~ 2 , values common to giant HII 
regions. 

We derived oxygen, sulphur, nitrogen, neon, and argon 
chemical abundances in knots B and C, where the auroral line of 
[Om] at 4363 A had been measured and thus the electron temper- 
ature of the high-excitation zone also was derived. We assumed 
also an inner structure of the electron temperature to deduce the 
values of the temperatures of [On], which is representative of 
the low-excitation zone, affecting the calculation of the chemi- 
cal abundances of + , N + , and S + . This inner ionization structure 
leads to derivation of the temperature of [Sra], valid for the inter- 
mediate zone (Garnett, 1992) and calculation of the ions of S 2+ 
and Ar 2+ . The derived ionic abundances in these two knots are 
very similar, leading to total O/H abundances of 12+log(0/H) = 
7.90 in knot C and 7.95 in knot B. These values deviate slightly 
from the metallicity reported in the literature for the integrated 
object, which is 12+log(0/H) = 8.24 (Shi et al., 2005; Kewley et 
al., 2005) and the difference is even greater than the value for the 
nuclear region, which could be considered as knot C, given by 
Kewley et al. (2005), which is 12+log(Q/H) = 8.55. The value 



reported by Shi et al. (2005) are based on the direct method, 
while those reported by Kewley et al. (2005) were derived using 
empirical calibrations since no temperature sensitive lines were 
observed. 

This is also the case in our observations of knots A and 
D. The calculation of metallicities in these two knots can only 
be carried out with calibrations based on the strongest emis- 
sion lines because no auroral lines were detected and the direct 
method cannot be used. 

The different strong-line methods for abundance derivations, 
which have been widely studied in the literature, are based on the 
directly calibrating of the relative intensity of some bright emis- 
sion lines against the abundance of some relevant ions present 
in the nebula. For the case of oxygen, we take the calibrations 
studied by Perez-Montero & Diaz (2005), who obtain different 
uncertainties for each parameter in a sample of ionized gaseous 
nebulae with accurate determinations of chemical abundances in 
the whole range of metallicity. 

In Figure [7] we show the corresponding total abundances as 
derived from several strong-line methods and the oxygen abun- 
dances calculated from the electron temperatures measured in 
knots B and C. Among the available strong-line parameters we 
studied the O23 parameter (also known as R23 and originally de- 
fined by Pagel et al. (1979) and based on [On] and [Ora] strong 
emission lines). This parameter is characterised by its double- 
valued relation with metallicity, with a very large dispersion in 
the turnover region. According to the values measured in knots 
B and C, we used the McGaugh (1991) calibration for the lower 
branch, obtaining similar values for the oxygen abundance in the 
four observed knots. 

The N2 parameter (defined by Storchi-Bergmann et al., 
1994) is based on the strong emission lines of [Nn]. It remains 
single-valued up to high metallicities in its relation to oxygen 
abundance, and it is almost independent of reddening and flux 
calibrations. Nevertheless, it has the high dispersion associated 
to the functional parameters of the nebula (ionization parameter 
and ionizing radiation temperature) and to N/O variations. We 
used the empirical calibration of this parameter from Denicolo 
et al. (2002) to derive the oxygen abundance in the four knots 
of this galaxy. We can see in Figure [7] that the abundances pre- 
dicted by this parameter are quite similar for all the knots al- 
though more than the values derived from the direct method in 
knots B and C. 

The S23 parameter was defined by Vflchez & Esteban (1996) 
and is based on the strong emission lines of [Sn] and [Sin] . The 
calibration done by Perez-Montero & Diaz (2005) yields oxy- 
gen abundances comparable to the four observed knots, slightly 
higher than the directly derived abundances of knots B and C, 
but still consistent with them within the errors. In the case of 
sulphur, both the directly and empirically derived abundances of 
knots B and C are slightly different with knot C showing a lower 
abundance. However it should be remembered that the near IR 
sulphur lines relative to the hydrogen recombination lines are 
more affected by reddening than in the case of oxygen when no 
Paschen lines are observed, which is the case. 

The parameter 03N2, defined by Alloin et al. (1979) de- 
pends on strong emission lines of [Om] and [Nn]. We used the 
calibration due to Pettini & Pagel (2004) and, as we can see in 
Figure [7] it has a very similar behaviour to that of N2. 

The combination of the S23 and O23 parameters gives 
S23/O23, defined by Diaz & Perez-Montero (2000), a parame- 
ter that increases monotonically with the oxygen abundance up 
to oversolar regime and which is very useful to study variations 
over wide ranges of metallicity (e.g. disks). We applied the cali- 
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Fig. 7. The oxygen abundances and their uncertainties for each observed knot, as derived using different empirical calibrators. From 
upper to lower panel: O23, N2, S23 and 03N2, at left and AT3O3 and S23/O23 at right. The lower panels of right column represent the 
sulphur abundance as derived from the S23 parameter and the N/O ratio as derived from the N2O2 parameter. The dash-dotted line 
represent the abundances and their uncertainties as derived from the direct method in knots B and C. 



bration found in Perez-Montero & Diaz (2005) and, in this case, 
there is almost a perfect coincidence with the values found by 
the direct method in knots B and C. Again, we find a similar 
value in knots D and A. 

The Ar3C>3 parameter, defined and calibrated by Stasihska 
(2006) as the ratio of [Arm] 7136 A and [Oin] 5007 A emission 
lines, predicts identical values for the metallicity in knots B and 



C, although noticeably higher than the directly derived ones. A 
rather high value is found in knot A and no data on the [Ann] 
line in knot D were available. 

Finally, the N2O2 parameter, defined by Perez-Montero & 
Diaz (2005) as the ratio between [Nn] and [On] emission lines, 
can be used to obtain the N/O ratio. Using this parameter, very 
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Fig. 8. Histogram of the distribution in mass fraction (left) and visual light (right) of the most probable stellar population models 
fitted by STARLIGHT for each of the observed knots. 



Table 6. Values of the extinction, total stellar mass, and fraction 
of the mass in stars younger than 10 Myr for each knot in the 
best-fit model predicted by STARLIGHT. 



ID 




log (M/M ) 


f(age < 10 Myr) 


A 


0.19 


6.13 ± 0.05 


0.08 


B 


0.46 


6.92 ± 0.05 


0.11 


C 


0.04 


7.52 ± 0.04 


0.008 


D 


0.25 


6.46 ± 0.04 


0.013 



little difference is found in N/O among the four observed knots, 
which are also consistent with the directly derived ratio. 

3.5. The stellar population 

The study of the stellar content in the four observed knots was 
carried out using STARLIGHT code, which calculates the com- 
bination of stellar libraries and the extinction law that repro- 
duce the spectral distribution of energy, to derive the properties 
of the stellar population in each of the knots better. The used 
Population Synthesis stellar libraries and the method to fit the 
observed spectra of each of the knots in the polar ring is de- 
scribed in Sect. 3.1. 

In Figure |8] we show the age distributions of the mass frac- 
tion (left) and the visual light fraction (right) for the four knots. 
All of them present a very young stellar populations with ages 
below 10 Myr, responsible of the ionization of the surrounding 
gas, although most of the stellar mass belongs to populations 
older than 100 Myr in almost all cases. The estimated total stel- 
lar mass and the fraction by mass of the stellar population with 
an age younger than 10 million years, responsible for the ioniza- 
tion of the gas, are listed in Table [6] for all the knots, together 
with the internal extinction, represented by the reddening con- 
stants, estimated by the model. The masses were corrected using 
an aperture factor, which takes into account the light not col- 
lected by the slit and which was calculated as the weighted mean 



of the ratios of the light in the whole region as measured in the 
Ha image and in the regions covered by the slit. 

Although it is not expected to find the same values of extinc- 
tion in the gas and the stellar population, there is a very good 
agreement between the values of the extinction estimated by 
STARLIGHT in knots C and D in relation to the values obtained 
from the decrement of Balmer. In knots A and B, the values are 
higher than those derived for the ionized gas. In these two knots 
the specific weight of the youngest population is higher. Since 
STARLIGHT assumes a common source of extinction for all the 
stellar populations entering the fit, the degeneracy between age 
and extinction tends to overestimate the extinction derived for 
the optimal fit. On the other hand, in knots C and D, which are 
probably contaminated by the light coming from the host galaxy, 
the contribution by the young stellar population is much lower, 
and the light at longer wavelengths is dominated by the older 
stellar population and the estimated average extinction is more 
representative. 



3.6. Ionising stellar populations 

Some properties of the emission knots can be obtained from 
the measured Ha flux, such as Ha luminosity, number of ion- 
ising photons, mass of ionising stars and mass of ionised hy- 
drogen (see Diaz et al. 2000). In order to obtain these quanti- 
ties, we analysed the Ha image of the galaxy, retrieved from 
the Palomar/Las Campanas atlas of blue compact dwarf galaxies 
(Gil de Paz et al., 2003). We defined elliptical apertures on this 
image for each of the four knots extracted in the spectroscopic 
observations and we measured all the flux inside the elliptical 
apertures up to the isophote corresponding to the 30% of the 
light in the brightest knot. These regions are shown in the right 
panel of Figure 1, and their sizes are given in the first column of 
Table [7] as the radius of a circular aperture of area equal to that 
encompassed by the refereed isophote. The observed Ha flux 
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Table 7. Derived properties of the observed knots using the extinction-corrected Ha fluxes measured in the images from Gil de Paz 
etal. (2003). 



ID 


Re 9 


I(Ha) 


L(Ha) 


Q(H°) 


M,,,„ 


M(HII) 


SFR 




(arcsec) 


10~ 14 erg • s _1 - cm -2 ) 


(10 38 erg ■ s- 1 ) 


(lO 51 ^- 1 ) 


(10 5 -M Q ) 


(10 4 -M Q ) 


M Q /yr 


A 


1.64 


1.25 


4.91 


3.59 


7.16 


3.64 


0.004 


B 


2.94 


4.83 


19.0 


13.9 


27.6 


14.1 


0.015 


C 


2.65 


3.96 


15.6 


11.4 


22.7 


11.5 


0.012 


D 


1.43 


1.10 


4.33 


3.16 


6.31 


3.21 


0.003 



was corrected in each knot for reddening using the values of the 
reddening constants, c(HyS), given in Table|3] The derived values 
are given in Table [7] 

The derived values of the masses of the ionising clusters 
can be compared with those provided by the STARLIGHT fit. 
STARLIGHT gives lower values than the ones we derive by fac- 
tors between 6 and 8 for knots A, B and, C. In the case of knot 
D, the difference is more than an order of magnitude. If we take 
the masses derived from the Ha fluxes, which are, in principle, 
lower limits since neither dust absorption nor a possible escape 
of photons are taken into account, and we use the calculated 
proportions of young to total mass given in Table [6] we obtain 
total masses between 0.9 x 10 7 and 4.9x 10 7 M for knots A, 
B and, D and 2.8x10 s M Q for knot C. This high value points 
to a large contribution to the older population by the bulge of 
the host galaxy. This is consistent with the continuum light dis- 
tribution shown in Figure Q] which is wider than that associated 
with cluster C and peaks at a region displaced from that of the 
cluster by 1.6 arcsec (about 144 pc). On the other hand for knot 
B the cluster shows a significant continuum of the same spatial 
extent as the line emission and the maxima of both distributions 
coincide. Knots A and D show very little continua of their own, 
so their underlying population probably corresponds to the outer 
parts of the host galaxy. 

Since the differences in metallicity among the different knots 
are not significant and the age distributions are also similar, this 
could indicate a common chemical evolution of these knots, 
probably related to the process of interaction with the compan- 
ion galaxy IIZw70. 

The star formation rate (SFR) for each knot was derived 
from the Ha luminosity using the expression given by Kennicutt 
(1998): 

SFR = 7.9 x 10~ 42 x L(Ha) 

. The derived values are also given in Table [7] They vary from 
0.003 M yr 1 for knot D to 0.015 M yr" 1 for knot B and pro- 
vide a total SFR for the ring of 0.035 M yr _1 , about half of the 
values quoted by Kewley et al. (2005) for the integrated object, 
which is 0.066 M yr _1 . Using the sizes of the emitting regions 
obtained from the Ha images and listed in Table [7j we obtain 
a rather constant SFR per unit area, of the order of 7xl0~ 8 M 
yr 1 pc 2 . This is higher than the average value of individual HII 
regions in polar rings given in Reshetnikov & Combes (1994) 
which is 3.2 xlO" 9 ™ 



' M yr" 1 pc 



-2 



3. 7. Kinematics and dynamics of the polar ring 

We have analysed the differential radial velocity along the slit us- 
ing the Ha and [Sn] 6717 A emission lines. The heliocentric ve- 
locities are displayed in Figure [9] superimposed on the observed 
Ha profile. We find an asymmetric rotation, as it was already 
stated by Reshetnikov & Combes (1994), possibly affected by 




spatial distribution (arcsecs) 



Fig. 9. Rotation curve of IIZw71. Filled circles represent helio- 
centric velocity derived from Ha and open squares from [Sn] 
6717 A. Solid line represents the spatial distribution of flux of 
Ha to be compared with the rotation curve. 



the expanding velocities of the bubbles of ionised gas surround- 
ing the knots of star formation which reaches 60 km/s in the case 
of knot C. 

Assumming that we see the polar ring of the galaxy edge- 
on, with a radial velocity of 85 km/s, and considering an optical 
radius of 20 arcsec, which is where we can measure the emis- 
sion lines with a good enough signal-to-noise ratio, we calcu- 
lated a dynamical mass inside this radius of (2.8 + 0.2) xlO 9 
M . We then calculated the L^ luminosity inside the radius of 
20 arcsec, internal to the polar ring, from the B brightness dis- 
tribution given in Cairos et al. (2001), obtaining L^ = 7.2xl0 8 
L . Therefore, we obtain a value of 3.9 for the M/Lb ratio inside 
the star forming-ring. This is close to the value of 2.8 found by 
Reshetnikov & Combes (1994) from optical observations within 
a distance of 30 arcsec from the centre. Considering that their 
considered blue light probably encompasses the light coming 
from the emission knots in the ring, that luminosity probably 
constitutes an upper limit and provides a lower limit to the Mj Lb- 



4. Summary and conclusions 

Spectrophotometric observations of the galaxy IIZw7 1 were car- 
ried out to study the physical properties of the bursts of star for- 
mation along its polar ring. We extracted information on the four 
brightest knots in Ha that we labelled from A to D. 

The electron temperature of [Om] was measured in knots B 
and C, allowing the direct derivation of ionic abundances of oxy- 
gen, sulphur, nitrogen, neon and argon. The total abundances of 
these species are in the same range of metallicities measured 
in HII galaxies, but they are slightly lower than the abundances 
previously reported from measurements of the integrated galaxy 
acording to strong-line calibrations. The metallicity in the other 
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two fainter knots, where the temperature sensitive line could not 
be detected, was estimated by means of different strong-line pa- 
rameters. In all cases, the estimated abundances are consistent 
with those derived for knots B and C by the direct method with 
the parameter involving the sulphur lines providing the abun- 
dance closer to it. The rest of the parameters slightly overes- 
timate the oxygen abundance. The N/O abundance, as derived 
from the N2O2 parameter (the ratio of the [Nn] and [On] inten- 
sities), is remarkably constant over the ring indicating that local 
pollution processes are absent. 

Although the underlying stellar population of the host galaxy 
is detected in the spectrum of the central knot (knot C), the simi- 
larity of the star formation histories in the four knots, as deduced 
from the fitting of SSPs, as well as the derived metallicities, point 
to a common chemical evolution of the polar ring. 

We calculated the SFR in the different knots studied from 
the Ha luminosities. The combined SFR for the ring amounts 
to half of the integrated SFR for this galaxy reported by previ- 
ous authors. The average SFR is also higher than the average 
reported value for HII regions in polar rings by more than an 
order of magnitude. 

Finally, for the differential velocity, we used the wavelength 
position of Ha and [Sn] to measure an assymmetric rotation of 
the ring with a mean value of 85 km/s at an optical radius of 
20". This gives a dynamical mass of (2.8 + 0.2) • 10 9 M and 
an M/Lb ratio of 3.9, close to the value reported previously by 
Reshetnikov & Combes (1994). The kinematics of the ring is 
significantly affected by the expanding bubbles of ionised gas, 
which in the case of knot C reaches 60 km/s. 
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